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The finite-density and finite-temperature phase portraits of the baby-Skyrme and Skyrme models 
are investigated. Both grand-canonical and canonical approaches are employed. The grand-canonical 
approach can be used to find the "natural" crystal structure of the baby-Skyrme model and it is 
shown to have triangular symmetry. The phase portraits include solid, liquid and phase-coexistence 
between solids and vacuum states. Furthermore, a chiral phase transition can be observed for 
both models. The phase portrait of the Skyrme model is compared with that expected of strongly 
interacting matter, and we also contrast our results to other models. 
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I. INTRODUCTION 

One of the great questions in the theory of strong inter- 
actions is the behavior of baryonic matter at finite den- 
sity and finite temperatures. Even though it is widely 
believed that quantum chromo-dynamics (QCD) is the 
correct theory to describe the strong interactions, its 
non-perturbative nature makes it extremely difficult to 
describe the low energy consequences of the theory. Lat- 
tice QCD 1^ is the preferred numerical technique to study 
the QCD Lagrangian directly. Here, QCD is modeled 
non-perturbatively by discretising the theory on a lat- 
tice, and a numerical approximation is made to the path 
integral describing either the Euclidean time evolution 
or the thermodynamics of the theory. The masses of 
mesons, baryons, and glueballs have already been cal- 
culated to reasonable accuracy 0. The limitations of 
the approach are both numerical (momentum cut-off, fi- 
nite lattice size), and more theoretical. The most im- 
portant difficulty is the description of fermions on the 
lattice. Since lattice QCD relies on using Monte Carlo 
methods, computing power is another of the key limi- 
tations, ff sufficient powerful computers were available, 
the inaccuracies due to discretising QCD could be re- 
duced by adopting larger gridsizes. Therefore, as com- 
puting power increases in future, more accurate results 
will be achieved. Even then, dealing with a finite den- 
sity of fermions is an especially tricky problem, since the 
fermionic determinant to be calculated becomes complex, 
leading to various instabilities. 

Due to these difficulties, it is quite common to use 
models to study the phase diagram of strongly interact- 
ing matter. Recent examples of such approaches can for 
instance be found in studies by Rajagopal et ^ , who 
describe the phase diagram using a Nambu-Jona-Lasino 
model [|| and a random matrix model Q. The QCD 
phase diagram, obtained from a model with two mass- 
less quarks in Ref. A typical phase diagram shows 
a hadronic-matter phase, the quark-gluon plasma, and a 
color superconducting phase. Earth's surface, in general, 
has a low temperature (0.025 eV) and a chemical poten- 
tial much less than 1 GeV, and therefore we live in a (liq- 
uid) hadronic-matter phase. The most striking feature 



of the phase diagram is the presence of a critical point. 
This separates an area where a first-order phase tran- 
sition occurs, from one where a continuous path can be 
taken from the hadronic phase to the quark-gluon plasma 
phase. At the critical point itself, it is believed that there 
is a second-order phase transition. The position of the 
critical point is extremely difficult to model and has been 
estimated to have a temperature Te ~ 140— f 90MeV and 
a chemical potential of he ~ 200 - 800 MeV |, §. The 
color superconducting phase cannot be modeled by the 
Skyrme model because color-flavor locking is observed, 
and in order to model it, three colors need to be consid- 
ered rather than the Nc ~ oo approximation underlying 
the Skyrme model, which will be discussed later. The 
color superconducting phase will therefore not be dis- 
cussed further. 

The Skyrme model has a venerable history j|] in the 
non-perturbative description of nucleon structure and the 
low-energy behavior of baryonic matter, since it contains 
a good description of the long-wave length behavior of 
the dynamics of hadrons. As has been argued by 't Hooft 
and Witten |l^ ^ , this is closely related to the large- 
number-of-colors limit of QCD, in which baryons must 
emerge as solitons, in much the same way as happens in 
the Skyrme model. Alternatively, we can interpret the 
model in terms of chiral perturbation theory, in terms 
what is much like a gradient expansion in terms of the 
pion field. The model can be used to describe, with due 
care [|l2, 13 , systems of a few nucleons, and has also been 
applied to nuclear and quark matter. Within the stan- 
dard zero-temperature Skyrme model description there 
are signatures of chiral symmetry restoration at finite 
density, but in a rather special way, where a crystal of 
nucleons turns into a crystal of half nucleons at finite 
density, which is chirally symmetric only on average [p^ . 
The question of the finite-temperature has never been ad- 
dressed and would be of some interest. We shall also look 
at the two-dimensional Skyrme model, which has the ad- 
vantage of being easier to interpret, but also has physical 
relevance; a special form of the two-dimensional Skyrme 
model has recently been developed for use in quantum 
Hall systems ||l5|. This model is obtained as an effective 
theory when the excitations relative to the v = 1 ferro- 
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magnetic quantum Hall state are described in terms of 
(a gradient expansion in) the spin density, a field with 
pro perties analogous to the pion field in nuclear physics 
|16| . Apart from obvious changes due to the number of 
dimensions, the new approach differs from the histori- 
cal Skyrme model by having a different time-dependent 
term in the Lagrangian, and the appearance of a non- 
local interaction, where the topological charge density at 
different points interacts through the Coulomb force. In 
the limit of large Skyrmions this last term can be ap- 
proximated by a more traditional local "Skyrme term" , 
which is quartic in the fields, leading to the standard 
baby-Skyrme model with local interactions. Even with 
the non-local complication the model is on the whole 
remarkably similar to the nuclear Skyrme model. The 
Skyrme-field effective degrees of freedom describe the 
ground state of such systems, and probably also the low 
energy dynamics and thermodynamics, so that we can 
ask similar questions about the finite-density physics as 
for baryonic matter. 

In this paper we shall consider the local baby Skyrme 
model, which can be regarded as the 2D version of the 
traditional 3D nuclear one. We can also extrapolate to 
one dimension. The model obtained in that way is the 
sine Gordon model expressed in terms of an unit com- 
plex field. All of these models contain a topologically 
conserved charge, which in the three-dimensional case is 
identified with baryon number. This gives rise to topo- 
logical solitons, which in nuclear physics are identified as 
baryons. Since the behavior of nuclear matter at high 
densities is expected to reveal, probably by one or two 
phase transitions, the substructure of baryonic matter, it 
is interesting to study this phenomenon in Skyrme mod- 
els. At the same time it allows us to look at the phase di- 
agram of quantum Hall Skyrmionic systems, which may 
well be the easiest way to calculate parts of the phase 
diagram of the underlying electronic model. 

At finite density but zero temperature the Skyrme 
models and the sine Gordon model have a crystalline 
structure, which consist of regularly spaced solitons. In 
the one-dimensional case, where exact solutions exist 
this is known to be a zero temperature artifact, and 
we have a liquid at any finite temperature, no matter how 
small. In the nuclear Skyrme model we would like to find 
a fluid, to mimic the quantum liquid behavior expected in 
nuclear matter jisj . In order to appreciate the subtleties 
involved, one must understand that the Skyrme model, 
as a classical field theory, can be understood as a semi- 
classical (large action) limit of a quantum theory. Glearly 
the quantum fluctuations in the underlying theory could 
be large enough to wash out the crystalline structure, as 
happens for the sine Gordon model. In the special case 
of one space dimension it is also well known that thermal 
and quantum fluctuations play exactly the same role, and 
thermal fluctuations also break the crystalline state. For 
the Skyrme and baby Skyrme models it is not easy to ac- 
cess the quantum fluctuations, since the field theories are 
non-renormalisable, but we can access the thermal fluc- 



tuations (with care, because nonremormalisability plays 
a role in the thermodynamics as well!). Furthermore, we 
might well wish to study the physics of these systems at 
temperatures where thermal fluctuations dominate the 
physics. Therefore, we shall concentrate on the finite 
temperature phases of these classical field theories. 

The most direct approach to the problem is to per- 
form a Monte-Carlo (Metropolis algorithm based) study 
of the partition function of each of these models. As we 
shall argue below the least obvious aspect in such an ap- 
proach is how to deal with the topological conservation 
laws. Rather than immediately tackle the 3D model, we 
shall first concentrate on the (local) 2D model. This has 
the advantage that the visualization of, and also the un- 
derstanding gained by, the results is much more straight- 
forward. The extensions to quantum-Hall Skyrmions are 
under way and will be presented elsewhere ]l9| . We ex- 
pect the current for the QHE solitons to be very similar 
to the results reported here. 

This paper is organized as follows. In Sec. II we discuss 



the formalism succinctly. In the next section. Sec. Ill we 



discuss in detail the results obtained in our simulations 
(a small excerpt has already been published in Refs. pC| , 
pi| ) . We draw some conclusions, and give an outlook for 
possible future work in Sec. 



II. FORMALISM 



A. Lagrange density 



The Lagrange density for Skyrme models can be given 
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+T(a^'^fc 3'^ 0i 5.0^9^00 + -^(00 - 1) , (1) 

which as long at we take the space index /i = 0, 1, . . . ,d 
to describe d+l dimensional space, and the field compo- 
nents k,l = 1, d + 1 to describe a d -I- 1 dimensional 
unit vector, can describe the Sine-Gordon theory {d — 1), 
the baby Skyrme model {d = 2) and the nuclear Skyrme 
model {d — 3). The Lagrange density can be separated 
into a kinetic and a potential term 



c = r[0] - v[0] 



where 



1 



T[<f>] = -dt<f>adt(t>bMab[ 



(2) 



(3) 



Mab[4'] = Sab + di(f>cdi(j)cSab ~ di(j)adi(j)t , (4) 



+ ^{[dz(pa di(l)a]'^ - [di(l)a dj(l)a][di(l)t dj(j)t]){5) 
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For our later calculation we need we need the expres- 
sion for the energy density, T + V, which for the Skyrme 
model becomes 

£ = \dt(i^a dtcbb Mab{ct>{x)) + V(0(a;)) . (6) 

The two sets of constraints (implemented by the delta- 
functions in the integrals) are cfP' — \ — Q and (f> ■ (p = 
at each point in space, where the first constraint forces 
to have unit length and the second constraint follows by 
taking the time derivative of the first constraint. 

In 3D the topological charge density, which was identi- 
fied by Skyrme with the baryon-number density operator, 
is given by 

llTT 

= ■^<iet{(f>,d^cf),dycf),d^(t}) , (7) 

with simpler expressions in ID and 2D. The energy den- 
sity (|^) and baryon-number density (|^) are then inserted 
into the grand-canonical partition function (||) . 



B. Partition EYinction 

In order to study the thermodynamics of a theory, we 
must construct a partition function. We shall discuss 
both grand-canonical and canonical partition functions, 
which are also discussed in Refs. [2^]. The explicit ex- 
pressions for the thermodynamic partition functions will 
be given for the baby-Skyrme and the Skyrme models. 
The Metropolis principle is used to evaluate the parti- 
tion functions, see Refs. M for more details. 



1. The Grand- Canonical Partition Function 

The general thermodynamic partition function for a 
general unitary field theory on the lattice is given by 

» M'' 



which we shall refer to as the scaled energy, is 



p=i 



X exp(-/5(£; - ^B)) , 



(8) 



where an integration is to be performed over the field 
and field derivative at each lattice point. The quantities 
E and B are the total energy and baryon number of the 
system. 

If the field makes small vibrations (which must be har- 
monic) at each lattice point then each degree of freedom 
contributes ^fc^T to the free energy. There are M'^ lat- 
tice sites, each with d potential and d kinetic degrees 
of freedom. The total energy due to lattice vibrations 
is thus E'vib — The energy due only to sohtons. 



scaled 



E< 



total 



E, 



total 



dM"^ 



(9) 



It is clear that these vibrational contributions are un- 
physical and need to be removed, since in the contin- 
uum limit their contribution would become infinite. The 
sine-Gordon model is renormalisable, and by subtract- 
ing the energy due to the lattice vibrations, we obtain 
the contribution due entirely to the solitons in a similar 
way as was done analytically in Ref. The baby- 

Skyrme and Skyrme models are non-renormalisable, and 
we find that it is therefore not possible to remove all 
the lattice-dependency from the energy expression. How- 
ever, by removing the harmonic contribution, the lattice- 
dependency is reduced. We will not use the energy of a 
system to examine phase transitions, because of these 
lattice-dependencies. Instead, we use other, less lattice- 



dependent methods, which we will discuss in Sees. IIIB 



and III C 



It is useful to use density-like quantities rather than to- 
tal quantities, and therefore we often use quantities like 
the average potential-energy density V = -ftdip- and the 
average baryon density B = j^ji^ where h'^M'^ is the 
volume of the system. Results are even more compara- 
ble when measuring ^ (or ^), which gives the potential 
energy per soliton. By measuring density-like quanti- 
ties, we are measuring quantities which remain finite and 
meaningful in the limit of an infinite simulation volume. 

The chemical potential /i determines the particle den- 
sity in a particular system. We initially expected that a 
simulation of the grand-canonical partition function au- 
tomatically includes the probability of adding or remov- 
ing a soliton from the system. Unfortunately, because 
we are working with extended particles, a soliton can 
only be added or removed from the system if a number 
of field vectors on our discretised lattice change simulta- 
neously. The probability of this happening is extremely 
small, and has never been observed in our simulations. 
The only examples of the total particle number changing 
is if the numerics become unstable. In this case, neigh- 
boring vectors are pointing in quite different directions 
such that the derivatives, and therefore the energy and 
baryon densities, are not calculated correctly. The nu- 
merical breakdown, i.e. the unphysical loss of winding 
number, is dependent on the lattice spacing, because it 
occurs at lower densities and temperatures if the lattice 
spacing is increased. Such numerical breakdown must 
occur when we have less than the minimum number of 
lattice points that can describe a Skyrmion. 

To solve the problem that the number of particles in 
a given simulation is fixed by the topology, we use an 
open system, where particles can flow through the bound- 
aries freely and therefore we have a heatbath that acts 
as a source of particles. The method we use to imple- 
ment our open system will be discussed in Sees. IIIB, 



and III C , when we discuss the thermodynamics of our 



soliton models. 
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2. The Canonical Partition Function 

The canonical approach is usually applied to a closed 
system, i.e. where particles are trapped in a box, and 
the number of particles in this box remain constant. As 
already mentioned in the last section, the number of soli- 
tons is conserved by the topology of the system. There- 
fore, the canonical approach seems to be very suitable to 
investigate the thermodynamics of soliton systems. The 
average density of particles is therefore defined by the 
initial choice of the number of particles B in the volume 
V and it never changes throughout the simulation. We 
choose not to implement fixed boundaries, which have the 
vacuum value on the simulation edges, but instead we use 
periodic boundaries. A system with periodic boundaries 
can be interpreted as an approximation of an infinite sys- 
tem because the simulation box is effectively duplicated 
an infinite number of times. Periodic boundary condi- 
tions imply that if a Skyrmion flows through one bound- 
ary of the simulation, then it reappears through the oppo- 
site side. Periodic boundary conditions have two advan- 
tages over fixed boundaries; firstly that the simulation 
can be interpreted as an infinite system, and secondly 
that we avoid unphysical boundary effects occurring from 
Skyrmions interacting with the fixed boundaries. 

The general expression for the canonical partition func- 
tion is 




X exp{-l3E) , (10) 

where all the symbols are again defined as in Sec. |II B 1 . 
The harmonic contributions to the potential energy are 
removed in a similar manner as for the grand-canonical 
approach, 

nM" 

Kscalcd = Vtotal • (H) 

Zp 

In Ref. [2^ the zero-temperature minimal-energy solu- 
tions solutions are calculated using simulated annealing 



within the canonical approach. For such static solutions, 
a renormalisation is irrelevant. Furthermore, the kinetic 
energy is irrelevant since the total energy is equal to the 
potential energy. 

We now use the fact that cf) is an eigenvector of the 
mass matrix (|4|) to perform the integral over cf) analyti- 
cally. The partition function now simplifies to 

X exp(-/3(Vp -/iSp)) , 
(12) 

where Vp and Bp are the potential energy and baryon 
number at the lattice site p and represent the potential 
energy and baryon number in the surrounding unit cell. 
The overall constant factor JV~^ is irrelevant when ap- 
plying the Metropolis principle, and will be ignored from 
now on. Eq. ( p^ will therefore be used as the grand- 
canonical partition function. The lattice does not need 
to be cubic, but we have adopted the convention that all 
the sides have equal length in order to avoid the addi- 
tional notation required when using unequal lengths. 

The procedure discussed above can be applied to the 
canonical formalism, and as expected we obtain the ex- 
pression ( |l^ ) with fi = 0. 

The partition function ( p^ ) cannot be evaluated ana- 
lytically, and therefore a Metropolis algorithm is applied, 
see Refs. Q for more details. 

The only exception is the sine-Gordon model. The 
grand-canonical partition function is derived in Ref. . 
It can very easily be compared to our simulations, since it 
is obtained by discretising the field and then taking the 
continuum limit. Like many other field-theory models, 
the observables tend to infinity as the continuum limit 
is being reached. It thus illustrates clearly the role of 
renormalisation. 

The canonical partition function of the sine-Gordon 
model in the angle representation is given by 



Zb = y n ^^'^P *^^P (^-Pi^dtUp dtttp + ^d^ap d^Up + 1 - cos ap)h^ 
j WS Up eyi^(^[3{]^dxapdxap + l- cos ap)h^ , 



const 1 



(13) 



where a is the field, and originates from evaluat- lattice spacing h is required to convert the energy density 

V ^ 1 to the energy in a cell from the discrete lattice. The con- 

ing the time derivative analytically, see Sec. |IBl[ The ^^^^^ ^^^^^^ expression for Zb is irrelevant when 
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measuring thermodynamic quantities, and is therefore set 
equal to one. The reason every lattice point behaves like 
a classical harmonic oscillator is that the kinetic energy 
and potential energy are quadratic in time derivatives 
and spatial derivatives, respectively. The contribution to 
the fluctuations at individual lattice points from the po- 
tential 1 — cos(q;) is negligible. The energy at a single 
lattice point due to fluctuations in the field is given by 



^ X /i, or ^ 



and h is the lattice spacing. Therefore, 



where i is the lattice index 



„ M / 

Hexpf 


i 


X dai . 


■ ■ dau 


const X 


. M/2 



{a 



2h 



(14) 



Using P — the harmonic oscillator contribution 



to the pressure becomes P = j;gin{- 



Since the 



harmonic oscillator at every lattice point gives a finite 
contribution to the pressure P, it becomes infinite in the 
continuum limit. The constant is not evaluated explicitly 
here because its contribution to the energy which we are 
interested in is irrelevant. Furthermore, the lattice spac- 
ing h within the logarithm also vanishes. The energy be- 
coming infinite in the continuum limit can be overcome 
by removing the contribution from the harmonic oscil- 
lators. This way, all quantities are calculated without 
contributions from the harmonic oscillators, and contain 
only contributions that are not divergent in the thermo- 
dynamic limit. Therefore, the pressure due entirely to 
the solitons is given by 



P' = - 



A 



8/32 



(15) 



We find that the density and internal-energy density 
for a given /? and /i, which are obtained by using the stan- 
dard thermodynamic relationships [^5[ ^ p — 

respectively, are 



and..-(^) 




FIG. 1: The exact (top figure) an numerical (lower figure) 
results for the grand-canonical equation of state of the sine- 
Gordon model. 



III. RESULTS 
A. Sine-Gordon model 

A check of the calculations can be made by comparing 
the exact results as discussed in the previous section to a 
calculation using open boundary conditions. This means 
that we impose no boundary conditions at all on the finite 
simulation volume, but allow solitons to enter the system 
through those boundaries. Even though, in principle, one 
can create topological charge on a grid, in practice topo- 
logical charge is very well conserved. The only way to 
perform a grand-canonical simulation in practice is thus 
using open boundary conditions. We then compare the 
simulated result for the density of solitons to the analytic 
resuh, Eq. (^, in Fig. | as a function of /3 and fi. As 
we can see the results are very similar, with some statisti- 
cal noise in the simulated results. Since it is much more 
straightforward to perform large simulations for the sine- 
Gordon model than for the 2D and 3D Skyrme models, 
this shows the best quality of results that can be expected 
(apart from a limitation in simulation time, since we have 
looked at quite a dense set of grid for the results). 



and 



dA 



^^,4/32 



df3 



8/32 dfi 



(16) 



(17) 



B. Thermodynamics of the Baby-Skyrme Model 

1. The Grand- Canonical Approach 

The special feature of the grand-canonical approach is 
that we use open-boundary conditions. This means that 
the number of particles present in a particular simulation 
can change when particles enter or leave the system over 
the boundary. The lattice points at a boundary behave 
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FIG. 2: The fi, p, T phase portrait for the baby-Skyrme 
model. Top: The order in which the points in the plot are 
calculated is in the increasing T direction with constant p. 
Bottom: The calculation order is in the decreasing T direction 
with constant p. The plots differ because the initial field 
configuration at each point is the final configuration of the 
previously calculated point. The axes are labeled in Skyrme 
units. 



slightly differently than the internal ones. If a field vec- 
tor is sampled at a boundary, it is accepted or rejected 
depending only on the change in the integrand on the 
neighboring internal plaquettes (we use the word plaque- 
tte for a square enclosed by four nearest neighbor lattice 
points). The field vectors at the boundaries are therefore 
less restricted than those at other lattice points. 

The p, p, T phase diagram for the baby-Skyrme model 
shown in Fig. y. Each point is a separate simulation 
which has reached equilibrium for a given chemical po- 
tential p and a temperature T. To save computing time, 
the initial field configuration used to calculate the den- 
sity p at each point in the plot was the final configuration 
for the previously calculated temperature. The arrows in 
Fig. H show the order the simulations where calculated in. 
Two plots of the same phase portrait have been shown 
where the order of calculation has been reversed. The 
reason for showing this is because the computing time re- 
quired to reach equilibrium becomes large near the phase 
transition, and the simulation time was not long enough 
to reach equilibrium near the phase transition. By com- 
paring the two graphs, one can see that it is difficult to 
predict the chemical potential of the transition to within 
5 Skyrme-energy units per soliton. In Fig. ^, the net den- 
sity of baby-Skyrmions p is shown against the chemical 
potential ^ at a fixed temperature T, for two different 
temperatures. In each plot, the order in which p changes 
has been shown; one in the increasing p direction and one 
in the decreasing p direction. For low temperatures, the 
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FIG. 3: The hysteresis effect when increasing and then de- 
creasing the chemical potential p. Top: The solid lines, where 
the arrows mark the direction of the simulation, show the hys- 
teresis effect for a low temperature T — 0.0001. The dashed 
lines are the results imported from Fig. ^ where the tem- 
perature was changed for constant p. Bottom: For higher 
temperatures, the hysteresis effect is less visible. The axes 
are labeled in Skyrme units. 



hysteresis effect is more visible than for higher tempera- 
tures, because solitons do not move as quickly and hence 
it takes longer to reach equilibrium. The magnitude of 
the hysteresis effect depends on the computing time used, 
as equilibrium is reached extremely slowly and we did not 
have the time to let the system strictly reach it. 

Each baby-Skyrmion has to have an energy greater 
than the topological lower bound for it to exist in an 
open system simulation. Fig. ^ shows that the chemical 
potential, which defines how much energy is given to each 
soliton, must be greater than approximately 18 Skyrme- 
energy units per baby-Skyrmion, p > 18. Once the chem- 
ical potential is greater than the threshold, the density 
increases approximately linear with the chemical poten- 
tial, p (X p- /itliroshold- 

As long as p is large enough, which is the case for al- 
most all of the finite-/? region, the state of the system is 
crystalline. The open boundaries allow crystal structures 
to form without defects, and therefore we can determine 
the "natural" crystal structure which will be described 
in the next section. The liquid state and the phase co- 
existence between solid and vacuum states exist near the 
phase transition and are too difficult to model using the 
grand-canonical approach. These states will be examined 



using the canonical approach in Sec. IIIB 2. 

For high densities and low temperatures, the 
Skyrmions merge into a lattice where it is not possible 
to identify individual Skyrmions, see Fig. ^. Since there 
are no fixed boundaries in this simulation, there are no 
defects in the crystalline solutions either. The lattice is 
determined most easily from the structure of the "holes" 
with zero baryon-number density. There are two such 
minima present per baryon. In this crystal phase, the 
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FIG. 4: Baryon-number density plots of the natural crystal 
structure of the baby-Skyrme model with open-boundary con- 
ditions. The baryon-number density is averaged over 20000 
steps. The axes are labeled in Skyrme-length units. Lighter 
shades represent higher baryon-number densities. The param- 
eters for this simulation are /x = 30, T = 0.05, and the number 
of solitons in the simulation volume is about 13 {rho = xxx). 

average field {4>) is zero, and therefore a mean-field clii- 
ral symmetry exists. 

The triangular nature of the crystal is evidence that 
the type of crystal favored by the baby-Skyrme model is 
independent of the lattice used to discretise the model. 
Even though the lattice is square, the crystal which forms 
does not align itself to it. Since this simulation allows the 
crystal structure having the lowest energy per baryon to 
form without being effected by boundary conditions, we 
call it the "natural" crystal structure. The field config- 
uration within a unit cell was copied (since it satisfies 
periodic boundary conditions) and simulated annealing 
with periodic boundary conditions was used to accurately 
investigate the zero-temperature properties. 

The correct high-density structure pO| where 
Skyrmions merge into a lattice is shown in Fig. ^. The 
lattice is determined most easily from the structure of 
the "holes" with zero baryon-number density, since it is 
not possible to identify individual Skyrmions. There are 
two such minima present per baryon. In the low-density 
phase, the average field {(pk) is polarized in the direction 
of the sigma field, but for high densities, the average field 
is zero. Therefore, a mean-field chiral symmetry exists 
in the high-density phase and not in the low-density 
phase. The energy per baryon ratio when minimizing 
the energy of a field configuration having the form of the 
natural crystal is 18.27 Skyrme-energy units at a density 
of 0.024 particles per cubic Skyrme-length unit. If the 
density is lowered further, then the minimal-energy 
solution consist of small numbers of Skyrmions bound 
together, i.e. the multi-Skyrmions. In fact, this can be 
defined as the point where the phase transition between 
the high-density and low-density matter occurs. 
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FIG. 5: The high-density phase at zero temperature. Top: En- 
ergy per baryon against density. Middle: A baryon-number 
density plot of the lowest energy per baryon number attain- 
able, i.e. the "natural" crystal structure. Lighter shades rep- 
resent higher baryon-number density. Bottom: A baryon- 
number density plot at high density. The baryon number is 
more evenly distributed throughout the structure. The mid- 
dle and right plot are drawn to the same scale. 

2. The Canonical Approach 

As we have argued before, with periodic boundary con- 
ditions it seems impossible to generate additional topo- 
logical charge. Nevertheless these calculations are much 
preferable over the open boundary calculations, since 
the results suffer much less from artifacts due to those 
boundaries (but the results will now depend much more 
strongly on the symmetries of the initial seed to the sim- 
ulation). This is still the preferred way to calculate the 
phase diagram. In order to create the p, T, ^ phase 
diagram, the lattice-dependent terms must be removed 
from the internal energy. Unfortunately, the lattice- 
dependency cannot be removed entirely, because the 
baby-Skyrme model is non-renormalisable. 

The internal energy per unit area u is found by using 
the thermodynamic relationship 

where is the simulation volume and B denotes the con- 
stant particle number defined by the initial conditions. 
Therefore, 
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where h is the lattice spacing and -^{V) is the potential 
energy per unit area. 

One can attempt to remove some of the lattice- 
dependency by assuming that the fluctuations at lattice 
points behave like harmonic oscillators, and therefore as- 
suming that the contribution from the potential energy 
of the harmonic oscillator, can be removed, so 

Mscalcd = -^(y)-^ . 

Unfortunately Uscaied still depends strongly on the lat- 
tice spacing. Furthermore, for high temperatures, Uscaicd 
can become negative, making its interpretation difficult. 
If the number of solitons in the system is zero, Uscaicd 
still has a non-linear behavior that can become nega- 
tive, when we would expect it to be zero. This is clearly 
an indication of the non-renormalisability of the Skyrme 
model, and we can no longer use the simple technique 
(count the lattice degrees of freedom) used to renormalise 
the sine-Gordon model 

Instead of searching for a detailed theoretical explana- 
tion of the lattice-dependency, which would require the 
study of counter terms, we chose to compare simulations 
with a fixed number of solitons with simulations that use 
the same lattice parameters and are at the same temper- 
ature, but contain zero solitons. Thus, we use 

Usc.icd = ^{V)-^{Vo) , (20) 

where (V) is the potential energy of a simulation with 
solitons, and (Vb) is the potential energy of the same 
simulation without any solitons. Although there is still 
some lattice-dependency, it is not very strong. It is likely 
that it is not possible to completely remove the depen- 
dency on the lattice spacing because the model is non- 
renormalisable. Nonetheless, we use this method because 
the internal energies for different lattice spacing become 
quite comparable, and phenomena such as phase transi- 
tion points do not seem to depend on the lattice scale 
chosen. 

We first use a canonical approach to study the p, T, ^ 
phase diagram. The internal energy per soliton is plotted 
against the density p and the temperature T in Fig. |^. 
The internal energy is measured using Eq. (|2^). The 
density p is altered by changing the lattice spacing h 
while keeping the number of Skyrmions in the system at 
a constant value. Although Eq. (|2^) is used to remove 
some of the lattice-dependency, there is still a remaining 
contribution. 




FIG. 6: The p, T, ^ phase diagram. Some of the contribu- 
tion from finite-lattice effects is removed by using Eq. ( |2C| ) 
to find the internal-energy density. The density p is altered 
by changing the lattice spacing and keeping the same number 
of Skyrmions in the system. The axes are labeled in Skyrme 
units. (16 solitons on a 90 x 90 lattice.) 



At finite temperatures, vibrations at each lattice point 
are present, but the contribution to the energy depends 
on the number of lattice points per unit volume. Since 
we vary the lattice spacing to change the soliton density, 
the internal energy per soliton as shown in Fig. |6|may not 
be properly compared for different densities. This phase 
diagram will not be discussed further, because it cannot 
be used to interpret the states of matter throughout the 
p-T plane. In the next two sections, the methods that 
successfully determine the states of matter for the baby- 
Skyrme model throughout the p — T plane are discussed. 

Along the lines first proposed by Klebanov ||2^, we 
analyse the grid-averaged fields (</)}. The average of the 
pion fields (tt) is always zero, but the average sigma 
field (a) changes. The average sigma field {a) is plotted 
against the temperature T and the density p in Fig. 0. 

When few Skyrmions exist in a large volume, {a) > 0. 
This is because the vacuum, where ct = 1, contributes 
a significant amount to the average sigma field. For 
high densities, however, the baby-Skyrmions combine in 
such a way that the individual baby-Skyrmions cannot 
be identified. Furthermore, an approximately chirally- 
symmetric configuration is formed where (cr) w and 
the sigma field and the pion fields are interchangeable 
without change in energy. For the baby-Skyrme model, 
chiral symmetry is never exactly satisfied (in a time av- 
erage), because the pion-mass term, which is required 
to stabilize the solitons, violates this symmetry. As the 
density of a system is increased, (cr) = is approached 



FIG. 7; The chiral symmetry (magnetization) (a) as a func- 
tion of p and T. For fluids and solids, we find (a) ~ 0. For 
tlie phase coexistence between solid and vacuum, (a) > 0. T 
and p are given in Skyrme units and (cr) is dimensionless. The 
parameter (a) has been measured from a single field configura- 
tion and has not been averaged through time, and the results 
show large fiuctuations. All simulations were performed for 
16 solitons on a 90 X 90 lattice, and the density was varied by 
changing the lattice spacing. 



asymptotically. In Fig. 0, a phase transition in (cr) can be 
observed to occur between the densities 0.02 < p < 0.04 
for various temperatures. 

The fluctuations in (cr) are largest in the region 0.02 < 
p < 0.08, T > 0.5. In the next section, we identify a 
liquid and a solid phase in the chirally symmetric phase, 
and the liquid region is located where the fluctuations in 
(cr) are largest. The liquid and solid states possibly ex- 
tend into the broken chiral symmetry phase, depending 
on where one defines the chiral phase transition to be. 
These fluctuations are present because in Fig. ^ (cr) was 
calculated from a single field configuration and not aver- 
aged in time. Also, the (cr) phase portrait seems to be in- 
dependent of the lattice parameters, unlike the internal- 
energy density. When doubling the number of lattice 
point and keeping the same number of baby-Skyrmions 
in the same volume, the (cr) phase diagram is not signif- 
icantly different than the one shown in Fig. 0. 

By examining the pseudo-time-averaged baryon- 
density plots at various temperatures and densities, one 
can easily identify three different states that a particu- 
lar simulation may be in, namely in a solid, in a liquid, 
or in a phase-coexistence state. Furthermore, correlation 
functions may be used to draw the same conclusions. In 
this section, the number of solitons in a system is varied 
to create different densities, while the lattice size and the 
lattice spacing remain fixed. Although the results shown 
have the same lattice parameters, the same conclusions 
can be drawn in simulations with different lattice param- 
eters, as has been checked for a number of examples. 

An example of the baryon-number density distribu- 
tion for a solid is given in Fig. ^. On the left we show 
a pseudo-time-averaged baryon-number density plot and 
on the right a snapshot where thermal fluctuations are 
evident. The approximate triangular crystal structure 



FIG. 8: A simulation in the solid phase. Pseudo-time- 
averaged baryon-number density plot. Lighter shades rep- 
resent higher baryon-number density and the units shown are 
in Skyrme-length units. The density p = 0.04, which implies 
that each light-colored circle contains two units of topological 
charge. 




FIG. 9; A simulation in a liquid phase, as a pseudo-time- 
averaged baryon-density plot with an average taken with 2000 
pseudo-time steps. The shades from dark to light represents 
low to high baryon-number densities and the units shown are 
in Skyrme-length units. As in Fig. ^ p = 0.04, but the tem- 
perature is higher. 



can be observed although defects are present, which are 
a consequence of using periodic boundary conditions, in 
contrast to the absence of defects seen when using the 
grand-canonical approach. 

A pseudo-time-average of a liquid shows an almost con- 
stant baryon-number density, because structures move, 
see Fig. ^ The visible structure in a pseudo-time- 
averaged baryon-density plot is dependent on the number 
of iterations used to create it, in this case 2000 pseudo- 
time steps. 

In the phase-coexistence region, one observes series of 
multi-solitons that do not change with pseudo-time, i.e. 
a solid structure surrounded by large regions of vacuum. 
These can either be interpreted as a crystal of nuclei, or 
as a percolating network. An example is shown in Fig. |l^. 

The solid, liquid and phase-coexistence states can also 
be identified by examining correlation graphs. The cor- 
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Solid 



FIG. 10: Pseudo-time averaged baryon-number density plot 
in the phase-coexistence region. The shades from dark to 
hght represents low to high baryon-number densities. As in 
Figs. pM p — 0.04, but the temperature is intermediate be- 
tween the two. 



relation graphs are created by plotting the correlation 
between the baryon number at the lattice points, 
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X (5(r- |r,-rb|), (21) 



against r, which is the separation between lattice points 
a and b. The correlation function C at r is calculated 
by taking the average of the correlation with respect to 
all a and b separated by a distance r. To keep track of 
the value B^B^ for each separation distance, we bin the 
values, and therefore the distances are rounded to the 
nearest half-lattice spacing, i.e. to the nearest h/2. 

The correlation graphs can be used to identify the state 
of a system. For solids and phase coexistence (which is a 
solid surrounded by vacuum), correlation functions show 
peaks a distance of several Skyrme-length units away 
from the origin, see Fig. |ll]. For liquids, however, the 
correlation function levels to zero quickly. The non-zero 
component is the correlation of Skyrmions with them- 
selves, and is therefore non-zero for approximately the 
length of the radius of a baby-Skyrmion. Some struc- 
ture is visible to a distance of 10 Skyrme-length units 
because interactions with neighboring Skyrmions are still 
present. It is because of the visible long-range interaction 
that we interpret the state to be a liquid rather than a 
gaseous state. The advantage of using correlation func- 
tions rather than pseudo-time-averaged baryon-density 
plots to identify liquids is that they are not dependent 
on the number of iterations used to create them. 

The transition from solid to liquid as the temperature 
is raised at a constant density is shown in Fig. [ij. Both 
the pseudo-time-averaged baryon-density plots and the 
corresponding correlation functions are shown to high- 
light the fact that the state of the system can be identified 
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FIG. 11: The correlation functions for a solid, a state in phase 
coexistence, and a liquid, respectively, (a): For the solid, there 
are peaks a distance of several Skyrme-length units away from 
the origin because the crystal is repetitive, (b): For phase co- 
existence, there are also peaks a distance of several Skyrme- 
length units away from the origin because the baby-Skyrmions 
are bound together, but the peaks are weaker than for a simu- 
lation in the solid phase, and lack the fine structure, (c): For a 
liquid, the correlation function decreases rapidly to zero after 
a few Skyrme-length units. The non-zero component exists 
because of the correlation of Skyrmions with themselves, (r 
is in Skyrme-length units and C is dimensionless.) 



through either method. The simulations whose baryon 
densities are plotted on the left have all been identified as 
solids and those plotted on the right have been identified 
as liquids. The pseudo-time average for these plots was 
set to 2000 iterations, and therefore it can be seen that 
the baryons move around more for simulations at higher 
temperature. When using more iterations, the baryon 
density approaches the average everywhere, but because 
there is less structure, it defeats the purpose of these il- 
lustrations. It is not possible to identify the structure at 
T = 0.2 directly from the baryon-density plot, but the 
correlation function clearly shows that it is a solid. The 
apparent motion of the solitons are lattice vibrations. 

Fig. ^ shows the final phase portrait for the baby- 
Skyrme model which has been created by examining the 
pseudo-time-averaged baryon-density plots and correla- 
tion functions at various temperatures and densities. We 
identify solid, liquid, and phase-coexistence regions, by 
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FIG. 12: Pseudo-timo-avcragcd baryon-dcnsity plots and the corresponding correlation functions for various temperatures with 
density p = 0.028. The length units shown are in Skyrme-length units and the correlation units are dimensionless. The contour 
levels for all plots are identical where white illustrates the highest baxyon density. All simulations are performed on a 90 x 90 
size lattice with lattice spacing h = 0.4. The simulations whose baryon densities are plotted on the left have all been identified 
as solids and those plotted on the right have been identified as liquids. 
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FIG. 13: The phase diagram. All plotted points are results 
obtained from simulations performed on a 90 x 90 lattice with 
a lattice spacing of /i = 0.4 Skyrme-length units. The axes 
are labeled in Skyrme units. 

combining the correlation function measure with that of 
the time-averaged baryon density. 

Unfortunately, we are unable to determine the order 
of any of the phase transitions. It is unlikely that any 
of them are of first order, because we would expect to 
see discontinuities in the p, T, ^ and p, T, {a) phase 
diagrams. Since higher-order phase transitions can only 
be identified by discontinuities in the n-th order deriva- 
tive, we are unable to do so since , the results presented 
here contain statistical fiuctuations that are large enough 
such that such an analysis becomes too difficult. 



3. Conclusions 

The baby-Skyrme model shows a rich phase diagram. 
By examining the chiral symmetry, the time-averaged 
baryon-density plots, the snapshot baryon-density plots 
and the correlation functions, it was possible to identify 
the state of matter that exists for a given density and 
temperature. Although the model is non-renormalisable 
and the internal energies are dependent on the lattice pa- 
rameters, all the other measured quantities used seem to 
be unaffected. Therefore, we are confident that the map- 
ping of the phase portrait of the baby-Skyrme model is 
given to a reasonable accuracy. 

A similar diagram will be created for the 3D Skyrme 
model, and the phase diagrams are expected to differ 
because the stabilizing terms have different forms. In 
particular, we would like to know at which temperatures 
and densities the phase transitions occur and compare 
these with putative QCD phase diagram such as those in 
Ref. @. 

C. Thermodynamics of the Skyrme Model 

In the previous section, we have successfully investi- 
gated the phase portrait of the baby-Skyrme model. We 
observed solid states, liquid states, and phase coexis- 



tence between solid and vacuum states. The internal 
energy, even after removing some of the lattice contri- 
bution, is still dependent on the lattice spacing because 
the baby-Skyrme model is non-renormalisable. However, 
by examining the chiral symmetry breaking, correlation 
functions, and baryon-density plots, we were able to de- 
termine the state of the system with methods that do 
not depend on lattice spacing. We will repeat the same 
type of analysis for the Skyrme model. Since the Skyrme 
model is often argued to be an approximation to QCD, 
the Skyrme phase portrait will be compared to other 
model of QCD (specifically Ref. which was created 
using another model). The phase transition between nu- 
clear matter and "quark" matter is of particular current 
interest. 

Unfortunately, we were unable to apply the grand- 
canonical approach to the Skyrme model. If the initial 
condition is a vacuum throughout the simulation box, we 
do not observe solitons entering the system. Skyrmions 
are expected to enter the system if /j, > 65 Skyrme-energy 
units per soliton. Simulations for a number of choices of 
/i, T and h suffered from numerical breakdown, which 
was identified by neighboring field vectors pointing in al- 
most opposite directions, where all derivative approxima- 
tions underlying the lattice approach break down. Nu- 
merical breakdown seems to occur more easily for the 
Skyrme model than for the baby-Skyrme model, most 
likely because the model is more complex. In the cases 
where a numerical breakdown did not occur, the initial 
vacuum never changed. Either solitons could not flow 
over the boundaries, or we did not have the computa- 
tional time required to observe any solitons entering the 
system. 

The field configuration that gives the lowest energy per 
soliton that has yet been observed is given by Castillejo 
et al and the energy per soliton is 3.8% above the 
topological lower bound. It is not clear whether this 
simple-cubic lattice in half- Skyrmions is the natural crys- 
tal structure. The periodic boundary conditions used in 
that reference favor this configuration. In fact, since we 
see that the baby-Skyrme model has a triangular pattern, 
we would expect that a hexagonal close-packed structure 
is a candidate for the natural crystal structure. We were 
not able to create a canonical simulation where we had 
the correct number of particles in a simulation box that 
was shaped to favor a hexagonal close-packed crystal. 
This does not imply that it is not the natural crystal 
structure, but instead it means that guessing the correct 
environment without knowing exactly what this crys- 
tal looks like is extremely difficult. The natural crystal 
structure cannot be identified from the canonical simu- 
lations either, because large simulation volumes are re- 
quired, and we do not have the computational resources 
to do this for the Skyrme model. Because of the large 
simulation volume required to identify the natural crys- 
tal structure from a canonical simulation, we probably 
would not have discovered the triangular crystal for the 
baby-Skyrme model without using the grand-canonical 
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FIG. 14: The p, T, ^ phase diagram. Some of the contri- 
bution from finite-lattice effects is removed by using Eq. ( |2^ ) 
to find the internal-energy density. The density p is altered 
by changing the lattice spacing and keeping the same number 
of Skyrmions in the system. The axes are labeled in Skyrme 
units. (32 Skyrmions on an 80 x 80 x 80 lattice.) 

approach. 

1. The Canonical Approach 

The canonical approach for the nuclear-Skyrme model 
is implemented in a similar manner as for the baby- 
Skyrme model. Again, some of the lattice-dependency 
is removed from the internal energy. The internal energy 
of the solitons is given by 

Uscaicd = ;^(V^)-;^(^o) , (22) 

where (V) is the average potential energy of a simulation 
with finite density and inverse temperature [3. The aver- 
age potential energy (Vb) is calculated from a simulation 
with identical parameters but without any Skyrmions, 
i.e. it is the potential energy due entirely to vibrations 
on the lattice. 

The internal energy per Skyrmion is plotted against 
the density p and the temperature T in Fig. |lj. The 
density p is determined by altering the lattice spacing h 
in a simulation with 32 Skyrmions on an 80 x 80 x 80 
lattice. Since some lattice-dependency is still present in 
the internal energy of the system, the internal energies 
at different densities may not be comparable, as they 
have been determined with different lattice spacings. It 
is not possible to determine the states of matter from 
Fig. |l^, and therefore alternative methods will be used 
in the following sections. 

The average sigma field is plotted against the den- 
sity p and the temperature T in Fig. |l^. As already 
discussed in the previous section for the baby-Skyrme 
model, a chirally-symmetric phase exists at high densi- 
ties, (cr) — 0, and chiral symmetry is broken at low den- 
sities, (cr) > 0. Although it is difficult to determine the 
position of the phase transition between low-density and 




FIG. 15: The p, T, (a) phase diagram. The kink is a sig- 
nature of the phase transition between the low-density and 
high-density phases. The axes are labeled in Skyrme units. 
(32 Skyrmions on an 80 x 80 x 80 lattice with varying lattice 
spacing.) 




FIG. 16: A baryon-density contour plot, with a high baryon- 
density contour, of a simulation at a high density (p — 
0.244 Skyrme-density units) and a finite temperature (T = 
0.2 Skyrme-energy units). An irregular, but crystalline, struc- 
ture can be seen. The mesh spacing is 0.08 Skyrme-length 
units. (64 Skyrmions on an 80 x 80 x 80 lattice with volume 
6.4"^ cubic Skyrme-length units.) 



high-density matter precisely, we believe it occurs where 
a kink can be observed near p = 0.2 for all temperatures 
plotted in Fig. |l^. 

The pseudo-time-averaged baryon-density plots, with 
the help of correlation functions, can be used to iden- 
tify the different states of matter. The solid, liquid and 
phase-coexistence states which are observed can then be 
used to determine the phase portrait for the Skyrme 
model. Since the Skyrme model is an approximation of 
QCD, the phase portrait we create can be compared to 
QCD. 

Examples of a contour plot from a pseudo-time- 
averaged baryon-density distribution for a solid are 
shown in Fig. |l6| for a high baryon-density contour. The 
crystal structure is irregular because the simulation vol- 
ume of 6.4'^ cubic Skyrme-length units is too small and 
therefore defects are introduced. We expected to observe 
a hexagonal close-packed crystal of half- Skyrmions, but 
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FIG. 17: A baryon-density contour plot of a simulation at 
a low density (p — 0.122 Skyrme-density units) and a finite 
temperature (T — 0.2 Skyrme-energy units). Eight 5 = 4 
multi-Skyrmions can be seen with temperature fluctuations 
and they would look like the B — 4 multi-Skyrmions. The 
contour shown for each the B — 4 multi-Skyrmion has an 
approximate volume of 2.8 cubic Skyrme-length units. The 
mesh spacing is 0.08 Skyrme-length units. (32 Skyrmions on 
an 80 X 80 X 80 lattice with volume 6.4^ cubic Skyrme-length 
units.) 

unfortunately no regular crystal can be identified. In a 
low-density contour plot we observe that the low- 
density regions are all connected. This is in contrast to 
the disconnected circles of vacuum observed for the baby- 
Skyrme model. 

We have not shown an example of a contour plot 
of the pseudo-time-averaged baryon-density distribution 
for a liquid. These look rather misleading when only 
one contour is shown, as regions with slightly higher- or 
lower-than- average baryon densities look like Skyrmions. 
These structures which are due to small differences 
in the pseudo-time-averaged baryon-number distribution 
are actually artifacts of the simulation process. When 
examining the results, we actually look at a variety of 
contours to verify that we have a liquid, but this would 
be too cumbersome to illustrate. A liquid can be dis- 
tinguished from a solid most easily by observing that a 
liquid does not have connected low-density contours that 
look like tubes. 

In the phase coexistence region, we observe multi- 
Skyrmions with large regions of vacuum surrounding 
them. An example is shown by a snapshot of a contour 
of the baryon-number distribution in Fig. |l^. Here, eight 
B = A multi-Skyrmions can be seen. The B ~ A multi- 
Skyrmions are very tightly bound, like the alpha parti- 
cles they model, which is why multi-Skyrmions of lower 
winding number are not seen. The boundary conditions 
also favor the B = A multi-Skyrmions over other multi- 
Skyrmions. Unfortunately, we are not able to run large 
enough simulations to see if a mixture of other multi- 
Skyrmions can be created. The regular spacing of the 
B = A Skyrmions is most likely also a consequence of a 
small number of Skyrmions in a limited simulation vol- 
ume. 
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FIG. 18: The correlation functions for a solid, a state in phase 
coexistence and a liquid, respectively, (a): For a solid, the 
correlations show long-range order at large distances. The 
correlation function corresponds to the contour plots shown 
in Fig. ^ (b): For phase coexistence, the self-correlations 
of the B — A Skyrmions are clearly visible for r < 1. The 
correlation between the different B = 4 multi-Skyrmions can 
be seen for r > 2.5. (c): The function falls off faster than 
the function for the solid, but medium-range interactions are 
still present and therefore it is not identified as a gas. The 
correlation function corresponds to the contour plot shown in 
Fig. (r is in Skyrme-length units and C is dimensionless.) 



C(r) 




FIG. 19: The correlation functions for Castillejo et al crystal, 
which is identified as a solid, (r is in Skyrme-length units and 
C is dimensionless.) 



The same correlation function that is used for the 
baby-Skyrme model can be applied to the Skyrme model, 
and has proven useful in the analysis of the results. The 
correlation functions for the solid, liquid and phase coex- 
istence states are given in Fig. The correlation func- 
tions seem to fall off faster for the Skyrme model than 
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FIG. 20: The phase diagram of the Skyrme model including 
a pion-mass term. The density of nuclear matter (NM) is 
marked at 0.16 fm~^. Other models predict the chiral phase 
transition at approximately 0.73 fm~^ All plotted points 
are results obtained from simulations performed on a 80 x 
80 X 80 lattice with a lattice spacing oi h = 0.08 Skyrme- 
length units. 



transition also appears to separate the phase coexistence 
region with the liquid and solid regions for this model 
at p = 0.18 Skyrme- length units for the temperatures 
shown. Fig. shows a sketch of the phase diagram and 
the states of matter and whether there is chiral symme- 
try is labeled. The solid and liquid phases are both chi- 
rally symmetric and are identified with the quark-gluon 
plasma. The coexistence phase is identified with nu- 
clear matter. The phase transition observed from chiral 
symmetry breaking and the phase transition observed by 
looking at average baryon-density plots and correlation 
functions are not necessarily signatures of the same phase 
transition. If they are two different phase transitions, 
they are close enough together so that we could not re- 
solve the two. Two different phase transitions may exist, 
because it may be possible to create crystals that break 
chiral symmetry which are not in phase coexistence with 
vacuum, which would be more identifiable with nuclear 
matter. 
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FIG. 21: A sketch of the phase diagram of the Skyrme model. 
The state of matter and chiral symmetry (or symmetry break- 
ing) is marked. The corresponding phases of QCD are marked 
in parentheses. 



for the baby- Skyrme model, but for solids and phase co- 
existence, a long-range order is clearly visible. The cor- 
relations for the putative solid must be compared to that 
of the parametrised solution by Castillejo et al, which 
consist of a simple cubic lattice of half-Skyrmions is 
shown in Fig. |l^. As one can see the decay of the correla- 
tion functions is very similar, if slightly less pronounced 
for the highly- defective solid found in our simulations. 
For liquids, there is also a visible medium-range order. 
The correlation of the 3 = 4 multi-Skyrmions with them- 
selves and with other B — -i multi-Skyrmions are clearly 
visible (at short and long distances respectively). Over- 
all, the correlation functions for the Skyrme model show 
the same type of behavior as for the baby-Skyrme model. 

The final phase portrait for the Skyrme model is shown 
in Fig. |2^. The different states of matter were distin- 
guished from each other by examining the pseudo-time- 
averaged baryon-density plots and the correlation func- 
tions. Classifying the different states of matter in the 
vicinity where the different phases meet is more diffi- 
cult than for the baby-Skyrme model, because long-range 
correlations are not as distinct. This means that the as- 
signments made at several points are somewhat ambigu- 
ous. Nevertheless, the approximate location where dif- 
ferent phases meet has been identified. The chiral phase 



2. Conclusions and comparison with QCD 

The phase portrait of the Skyrme model including a 
pion-mass term has been determined to a reasonable ac- 
curacy. The difficulties were mainly caused because each 
simulation required a significant amount of computing 
time. Furthermore, analysing three-dimensional struc- 
tures also requires more time and effort. Unfortunately, 
we did not have the resources to extend the energy per 
baryon and chiral symmetry-breaking phase diagrams to 
higher temperatures. 

We will now compare our results from investigating the 
Skyrme model at finite temperatures and finite densities 
to the phase portraits of QCD obtained by other models 
and illustrated in Sec. I. The phase transition between 
nuclear matter and deconfined quark matter occurs at 
p — 0.18 Skyrme-density units (2.9 fm""^). This phase 
transition, which has the characteristic that chiral sym- 
metry is broken, has been observed at p = 0.045 Skyrme- 
density units (0.73 fm~^), which was obtained from 
Ref. 1^. The phase transition occurs at about 18 times 
the density of nuclear matter on earth of approximately 
0.0099 Skyrme-density units (0.16 fm'^) [|||. The low- 
est energy per baryon for QCD is most likely at the nor- 
mal density of nuclear matter, because any excess en- 
ergy would be radiated away. For the Skyrme model, 
the lowest energy per Skyrmion occurs at 0.043 Skyrme- 
density units (0.70 fm""^). The lowest energy per baryon 
for the Skyrme model and for QCD are not easily com- 
pared, because this value will change significantly when 
the Skyrme model is quantised. 

The lowest temperature where a liquid phase is ob- 
served is at 0.5 Skyrme-energy units (8.45 MeV). Al- 
though this is approximately the temperature of the crit- 
ical point where nuclear matter can continuously (i.e. 
without a phase transition) transform from a liquid to 
a gas, we believe that the lowest temperature where a 
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liquid can be observed should be identified with the tem- 
perature of the critical point where nuclear matter can 
transform to a quark-gluon plasma continuously. These 
critical points have been obtained from Ref . Q . The rea- 
son we identify the lowest temperature at which a liquid 
is observed with the temperature of the critical point is 
because the high-density Skyrme matter, which behaves 
like a solid, is identified with the quark-gluon plasma, 
as individual Skyrmions are no longer distinct. The liq- 
uid phase may be the state where high-density Skyrme 
matter (quark-gluon plasma) can continuously be trans- 
formed to low-density Skyrme matter (nuclear matter). 
If this is the case, then the predicted temperature of the 
critical point is different by a factor of ten! 

Alternatively, the liquid state for the Skyrme model 
may be identified with the liquid state for the quark- 
gluon plasma, because the liquid state observed for the 
Skyrme model is chirally symmetric. (Since there are not 
many samples for the nuclear- Skyrme model, we observe 
that this is also the case for the baby-Skyrme model.) If 
this is the case, we have not been able to heat simulations 
enough to reach the critical point. Furthermore, it is 
unclear whether the Skyrme model shows any observable 
phenomena there. Possibly there is a liquid state for all 
densities above this temperature. 

The location of the phase transitions observed for the 
Skyrme model will not correspond exactly with those of 
QCD, because the Skyrme model is only a semi-classical 
approximation of QCD. However, it is difficult to model 
the different phases of QCD by using only one model. In 
future, the Skyrme model may be used to make state- 
ments about the order of phase transitions in QCD. We 
were not able to investigate this because we were only 
able to generate few samples (which each required a con- 
siderable amount of computation time) and even these 
samples are not statistically accurate enough to take 
derivatives of the phase spaces that were generated. We 
have, however, shown that field theories containing soli- 
tons have an interesting phase portrait. 

IV. CONCLUSIONS AND OUTLOOK 

The phase portraits for the baby-Skyrme and nuclear- 
Skyrme models were created accurately using the canon- 
ical approach. We have shown that the low-density and 
high-density phases can be distinguished through chiral 
symmetry breaking. Furthermore, pseudo-time-averaged 
baryon-density plots and correlation functions can be 
used to identify the state of matter at a given temper- 
ature and density. The baby-Skyrme model has a rich 
phase diagram. Although the model does not have a 
physical interpretation, we believe that similar phases ex- 
ist for Skyrmions in the quantum Hall effect. The phase 
diagram for the nuclear-Skyrme model was compared to 



other models of the QCD phase diagram. We were able 
to identify the phase where chiral symmetry is broken 
with hadronic matter and the chirally symmetric phase 
with the quark-gluon plasma. The location of the phase 
transition is at a higher density than predicted by other 
models. It remains unclear whether the boundary of the 
phase coexistence region with the liquid and solid regions 
occurs on the same line as the chiral phase transition. 

The grand-canonical approach could not successfully 
be used to generate accurate phase portraits for the 
baby-Skyrme and nuclear-Skyrme models, but it is use- 
ful in determining the natural crystal structure for the 
baby-Skyrme model. The open boundaries in this ap- 
proach allow crystals to form without defects. The grand- 
canonical approach is already being applied to QHE sim- 
ulations p9| . Unfortunately, this method was unsuc- 
cessful for the nuclear-Skyrme model, partly because we 
did not have the necessary computer resources. We also 
were not able to run simulations for a large number of 
different parameters. We believe that it is possible to 
determine the natural crystal structure for the nuclear- 
Skyrme model, and therefore propose that this be at- 
tempted again. We suggest to remove the pion-mass 
term in such calculations, because the minimal energy 
per baryon occurs in the high-density phase rather than 
the low-density phase (as is the case when a pion-mass 
term is included). When using the grand-canonical ap- 
proach, it may be effective to use an approximation to 
the derivative that results in overestimating the baryon 
number, rather than underestimating like we have done. 

The thermodynamic investigations can also be applied 
for the Skyrme model with a sixth-order term and for 
the w-stabilized Skyrme model. The higher-order terms 
become dominant at high densities, and may therefore 
be better models for describing such phases. 

The nuclear-Skyrme model is a semi-classical approx- 
imation to QCD in the Nc — > oo limit. This approxi- 
mation is not particularly accurate, but it means that 
a single model can be used to describe matter at both 
low densities and high densities. We were not able to 
investigate the order of the phase transitions, and since 
this may be relevant to QCD, it would also be interesting 
subject for future research. 

Thus, we have derived methods to investigate models 
with soliton solutions for both zero and finite tempera- 
tures. The methods we developed can easily be applied 
to a wide range of similar models. 



Acknowledgments 

OS has been supported by an EPSRC studentship; 
NRW acknowledges support from the EPSRC through 
research grants GR/L22331 and GR/N15672. 



17 



[1] I. Montvay and G. Miinster, Quantum Fields on a Lattice 

(Cambridge University Press, Cambridge, UK, 1994). [16] 

[2] R. Gupta, Introduction to Lattice QCD (????), Lectures 

given at the LXVIII Les Houches summer school, "prob- [17] 

ing the stand ard model of part icle interactions" , July 28- 

Sept 5, 1997. tiep-lat/980702^ . [18] 

[3] K. Rajagopal, Nucl. Phys. A 661, 150 (1999). [19] 

[4] M. Stephanov, K. Rajagopal, and E. Shuryak, Phys. Rev. [20] 
Lett. 81, 4816 (1998). 

[5] J. Berges and K. Rajagopal, Nucl. Phys. B 538, 215 
(1999). 

[6] M. A. Halasz, A. D. Jackson, R. E. Shrock, M. A. [21] 
Stephanov, and J. J. M. Verbaarschot, Phys. Rev. D 58, 
096007 (1998). [22] 
[7] K. Rajagopal (2000), hep-ph/0005101v2. [23] 
[8] T. H. R. Skyrme, Nucl. Phys. 31, 550 (1962). [24] 
[9] G. 't Hooft, Nucl. Phys. B 72, 461 (1974). 
[10] G. 't Hooft, Nucl. Phys. B 75, 461 (1974). [25] 
[11] E. Witten, Nucl. Phys. B 223, 433 (1983). 

[12] E. Braaten and L. Carson, Phys. Rev. D 38, 3525 (1988). [26] 
[13] E. Braaten, S. Townsend, and L. Carson, Phys. Lett. B 

235, 147 (1990). [27] 
[14] A. D. Jackson and J. J. M. Verbaarschot, Nucl. Phys. A [28] 

484, 419 (1988). 

[15] S. M. Girvin, in Topological aspects of low dimensional [29] 
systems. Les Houches lectures, edited by A. Comtet, 
T. Jolicoeur, S. Ouvry, and F. David (Springer Verlag, 



Berlin, 1998), vol. 29, p. 53, cond-mat/98. 

D. H. Lee and C. L. Kane, Phys. Rev. Lett. 64, 1313 

(1990). 

N. Gupta and B. Sutherland, Phys. Rev. A 14, 1790 
(1976). 

T. D. Cohen, Nucl. Phys. A 495, 545 (1989). 
T. Weidig and N. R. Walet (2001), cond-mat/0106157. 
O. Schwindt and N. R. Walet, in Intersections of Parti- 
cle and Nuclear Physics, edited by Z. Parsa and W. J. 
Marciano (AIP Conference Proceedings, Vol. 549, New 
York, 2000). 

O. Schwindt and N. R. Walet, Europhys. Lett. 55(5), 633 
(2001). 

T. H. R. Skyrme, Nucl. Phys. 31, 556 (1962). 

0. Schwindt, Ph.D. thesis, UMIST (2001). 

M. Hale, O. Schwindt, and T. Weidig, Phys. Rev. E 62, 
4333 (2000). 

F. Mandl, Statistical Physics (John Wiley, New York, 
1988). 

M. H. Everdell, Statistical Mechanics and its Chemical 
Applications (Academic Press, London, 1975). 

1. Klebanov, Nucl. Phys. B 262, 133 (1985). 

L. Castellejo, P. S. J. Jones, A. S. Jackson, and J. J. M. 
Verbaarschot, Nucl. Phys. A 501, 801 (1989). 
T. S. Walhout, Nucl. Phys. A 484, 397 (1988). 



